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PREFACE 


This  quarterly  report  was  prepared  by  Hughes  Research  Labora- 
tories, Malibu,  California  under  Contract  F30602 -75 -C -0001 . It  describes 
w°rk  performed  from  1 January  1975  to  31  March  1975.  The  principal 
investigator  and  principal  scientist  is  Dr.  James  E.  Pearson.  The  project 
is  part  of  the  adaptive  optics  program  in  the  Optoelectronics  Department, 
managed  by  Dr.  Viktor  Evtuhov,  at  the  Hughes  Research  Laboratories. 


SUMMARY 


The  objectives  of  this  research  program  are  threefold:  (1)  to 

experimentally  study  the  ability  of  a multidither  COAT  system  to  compensate 
for  thermal  blooming  and  turbulence,  particularly  when  it  operates  against 
complex,  multiple -glint  targets;  (2)  to  use  computer  simulation  to  cross 
check  the  experimental  results  and  to  investigate  areas  beyond  the  capabilities 
of  the  experimental  hardware;  and  (3)  to  perform  preliminary  studies  of  the 
effects  of  speckle  noise  modulation  on  COAT  system  convergence  and  per- 
formance. This  report  summarizes  the  work  performed  during  the  third 
quarter  of  the  contract  from  1 January  1975  to  31  March  1975. 

We  have  experimentally  studied  the  effects  of  target  motion-induced 
speckle  modulation  on  the  D ARP  A/ R ADC  COAT  system.  A scaled  laboratory 
experiment  using  a rotating  scotchlite  target  produced  a wideband  speckle 
noise  spectrum  over  the  dither  frequency  band  which  did  not  observably  inter- 
fere with  the  convergence  properties  of  the  system.  In  another  series  of 
experiments,  an  acousto-optic  modulator  with  64%  modulation  efficiency  was 
placed  in  the  DARPA/RADC  COAT  optical  system  and  swept  over  a variable 
frequency  bandwidth  in  a time  comparable  to  the  servo  loop  convergence  time. 
The  speckle  modulation  artificially  generated  by  this  technique  produced 
interference  signals  over  the  dither  band  greater  than  the  normal  dither 
signals,  and  produced  converged  power  degradation  of  35%.  The  results 
obtained  thus  far  indicate  that  the  system  performance  can  be  significantly 
degraded  if  the  modulation  is  large  enough,  but  the  systematic  evaluation  of 
which  targets,  target  ranges,  and  target  motions,  if  any,  can  actually  produce 
such  large  modulations  will  be  pursued  on  another  program. 

Computer  simulation  modeling  of  COAT  system  performance  in  the 
presence  of  speckle  noise  modulation  has  been  used  to  obtain  the  average 
convergence  level  of  an  18-channel  multidither  COAT  system  in  the  presence 
of  this  noise.  An  oscillation  of  the  system  about  the  average  convergence 
level  has  been  observed.  A statistical  analysis  of  the  COAT  system  and  of 
the  speckle  noise  has  been  developed  for  use  in  explaining  the  simulation 

results. 
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Multidither  adaptive  tracking  and  focus  (Beam  Active  JTrack  or  BAT) 
controls  have  been  implemented  into  the  DARPA/RADC  COAT  system  and 
into  a computer  simulation  code.  The  experimental  tracking  controls  have 
demonstrated  a convergence  time  of  8 to  10  msec,  the  design  goal.  A 
tendency  for  this  system  to  lock  onto  sidelobes  of  the  transmitter  array  has 
been  observed.  This  behavior  can  be  explained  in  terms  of  the  sidelobes  of 
the  18-element  COAT  planar  array  and  in  terms  of  the  tracking  dither  fre- 
quency being  larger  than  the  open-loop  bandwidth  of  the  phase-control  servo. 
The  BAT  computer  simulation  uses  2-axis  tracking  and  2-axis  focusing 
(2  separately  controlled  cylindrical  lenses).  For  moderate  turbulence 

(C^  Z = 4.  6 x 10~  1 1 m1  /3),  the  BAT  simulation  shows  an  improvement  in 
N 2 

average  strehl  ratio  from  0.67  to  0.86.  Fcr  strong  turbulence  (CNZ  = 21  x 

10  ^ ^) , the  improvement  is  from  0.2  to  0.5.  Almost  all  of  this 

improvement  can  be  produced  by  the  2-axis  focus  controls  alone;  adding  the 

tracking  controls  to  the  2-axis  focus  provides  only  a small  increase  in 

strehl  ratio.  A 37 -channel  multidither  COAT  system  achieves  strehl  ratios 

of  0.95  and  0.76  for  these  two  turbulence  levels. 

Experimental  observations  with  artificial  turbulence  equivalent  to 
C^Z  = 5 x 10'11  m1;  3 have  shown  good  agreement  with  the  simulation  data. 

We  have  observed  an  improvement  in  strehl  ratio  from  0.65  to  0.87  using 
only  the  tracking  controls.  When  the  18-channel  phase  controls  are  also 

used,  the  strehl  ratio  increases  to  0.95. 

A flowing  gas  cell  has  been  built  and  tested  for  thermal  blooming 
studies.  With  the  absorption  region  occurring  in  the  last  81%  of  the  focused 
propagation  path,  the  18  phase  controls  plus  two-axis  tracking  controls 
produced  a 35%  increase  in  peak  target  irradiance.  This  correction  is  still 
less  than  expected  and  further  study  is  indicated. 

An  exact  simulation  of  both  a multidither  servo  and  time -dependent 
thermal  blooming  has  been  determined  to  be  very  costly,  so  no  such  computer 
runs  have  been  made.  We  have  shown  that  finite  size  target  glints  together 
with  return-wave  COAT  systems  can  reduce,  rather  than  increase  target 
irradiance  in  the  presence  of  blooming.  These  effects  are  not  expected  to  be 
important  with  outgoing -wave  multidither  COAT  systems  (such  as  the 
DARPA/RADC  system). 
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During  the  final  contract  quarter,  work  will  be  concentrated  on 
experimental  measurements  of  simultaneous  thermal  blooming  and  turbulence. 
Analytical  work  will  focus  primarily  on  comparing  computer  simulation 
results  on  speckle -modulation  effects  to  experimental  observations  and  on 
using  the  statistical  models  we  have  developed  to  explain  the  observations. 


I. 


INTRODUCTION 


A.  Program  Objectives 

There  are  two  primary  objectives  of  this  program.  The  first  objec- 
tive is  to  determine  the  performance  limits  of  multidither  coherent  optical 
adaptive  techniques  (COAT)  by  performing  scaled  laboratory  experiments 
which  are  designed  to  produce  quantitative  data  on  the  nature  of  thermal 
blooming  and  turbulence  and  on  the  ability  of  multidither  COAT  to  correct 
for  these  distortions,  particularly  when  the  target  has  many  moving,  time- 
varying  glints.  The  second  objective  is  to  use  computer  simulation  as  a 
cross  check  on  the  experiments  and  as  an  analytical  tool  to  extend  the  under- 
standing of  propagation  distortions  and  of  COAT  systems  beyond  the  capa- 
bilities of  the  experimental  hardware. 

B.  Research  Program  Plan 

The  performance  of  multidither  COAT  has  already  been  proved  with 
an  18-element  visible  wavelength  system  developed  and  tested  on  DARPA/ 
RADC  Contract  F30602 -73 -C -0248  which  was  concluded  in  July  1974.  This 
same  system  is  used  for  the  laboratory  experiments  in  this  program.  The 
computer  simulation  codes  for  atmospheric  turbulence,  thermal  blooming, 
and  the  COAT  system  are  being  developed  on  DARPA/NOL  Contract  N60921- 
74-C-0249.  The  previous  DARPA/RADC  COAT  contract  also  supplied  a 
design  for  a gas  absorption  cell  which  has  been  used  to  simulate  convection- 
dominated  thermal  blooming  in  the  atmosphere. 

The  research  program  for  this  contract,  illustrated  in  Fig.  1,  runs 
from  2 July  1974  through  30  June  1975.  The  required  oral  presentation  was 
made  on  3 December  1974  as  part  of  the  DARPA/NOL  Adaptive  Optics 
Symposium  held  at  Lincoln  Laboratories.  ^ A contract  amendment  negotiated 
during  the  second  contract  quarter  provided  for  the  addition  of  auto -tracking 
and  autofocus  controls  to  both  the  DARPA/RADC  COAT  hardware  and  the 
Naval  Ordnance  Laboratories  computer  code.  A second  amendment,  nego- 
tiated during  the  third  quarter,  eliminated  a high-power  design  guidelines 
task  so  that  time  and  funds  could  be  directed  toward  experiments  on  COAT/ 
target-signature  interactive  effects.  The  schedule  in  Fig.  1 reflects  these 
changes. 
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II. 


TECHNICAL  ACCOMPLISHMENTS 


A.  Adaptive  Tracking  and  Focus  Controls 

Contract  Amendment  No.  1 provides  for  supplemental  funding  from 
the  Air  Force  Weapons  Laboratory  to  add  multidither  tracking  and  focus 
controls  both  to  the  DARPA/R ADC  experimental  COAT  hardware  and  to  the 
Hughes  multidither  COAT  computer  simulation.  The  controls  will  be  used 
both  separately  and  in  conjunction  with  the  DARPA/RADC  18-channel  phase 
controls  to  study  turbulence  and  thermal  blooming  compensation. 

In  the  following  discussions,  we  will  refer  to  three  possible  COAT 
system  configurations.  The  term  "18-channel"  will  be  applied  to  the  18- 
element  planar  array  which  produces  piston-type  phase  control  on  18  separ- 
ate transmitted  beam  segments.  The  acronym  "BAT"  (defined  below)  will 
apply  to  the  new  tracking  and  focus  controls  alone.  The  term  "21-channel" 
will  be  used  whenever  the  BAT  controls  are  used  together  with  the  18  phase 
control  channels. 

1.  Experimental  Results 

A multidither  tracking  and  focus  control  system  has  at  various  times 
been  referred  to  by  many  names:  "Conical  Scan"  (tracking  only),  "Beam 

Active  Tracking"  (BAT),  or  "Adaptive  Laser  Optical  Techniques"  (ALOT). 
For  simplicity,  we  will  refer  to  the  controls  developed  here  as  BAT.  The 
three -channel  BAT  system  built  for  the  DARPA/RADC  COAT  system  is 
shown  in  Fig.  2.  The  tracking  controls  are  designed  to  have  a 50  Hz  open- 
loop  unity  gain  bandwidth  and  a 10  msec  convergence  time.  They  use  two 
dither  frequencies  at  750  Hz  and  1000  Hz  or  a single  dither  at  1000  Hz,  with 
the  sine  signal  applied  to  one  channel  and  the  cosine  to  the  other  (conical 
scan).  The  focus  channel  is  a factor  of  10  slower  using  a 100  Hz  dither  and 
designed  to  converge  in  100  msec  (5  Hz  unity  gain  bandwidth). 

The  focus  control  shown  in  Fig.  i uses  a single  piezoelectric  (PZT) 
cylinder  to  deform  a single  mirror  into  a spherical  shape.  The  unit  was 
designed  to  use  a 2.25  in.  long  PZT5-H  cylinder  and  to  have  a loaded  deflec- 
tion  sensitivity  of  250  V/pm.  We  have  been  unable  to  obtain  PZT5-H  in 
the  desired  dimensions  and  so  have  used  PZT-4.  With  this  material,  a 
deflection  sensitivity  of  520  V/pm  is  expected. 


/ 


V 


1 ? 


MIRROR 


METAL  TAB 


SPACER 


ELECTRICAL 

GROUND 

CLIPS 


P2T 

CYLINDER 


PRELOAD 
END  PLATE 


HIGH  VOLTAGE 
LEAD 


Static  deflection  tests  of  the  focus  actuator  in  a Twyman-Green 
interferometer  are  shown  in  Fig.  4.  The  observed  deflection  is  quite  uniform 
and  spherical.  The  observed  sensitivity  is  1260  V/pm,  a factor  of  2.4  lower 
lhan  the  design  value.  We  are  not  sure  why  the  observed  sensitivity  is  so 
much  lower  than  expected.  One  possibility  is  that  the  mirror  spring  constant 
is  larger  than  our  calculations  pro  licted  relative  to  the  PZT  spring  constant. 
The  initial  curvature  (V  = 0)  is  an  adjustable  preload  on  the  PZT  cylinder  to 
provide  a spring  for  the  actuator  to  work  against  and  to  eliminate  the  need 
for  a bond  between  the  cylinder  and  the  mirror.  The  lowest  resonant  fre- 
quency of  this  device  was  determined  to  be  3.4  kHz. 

Because  of  some  difficulties  in  the  COAT  system  optics,  we  have  not 
used  the  focus  actuator  in  the  system  yet.  During  the  next  quarter  we  expect 
to  resolve  these  difficulties  so  that  the  calibration  and  propagation  tests  with 
this  device  can  be  completed. 

The  tracking  controls  have  been  tested  with  the  COAT  system  and 
some  preliminary  blooming  compensation  studies  are  reported  in  Section  II-C. 
The  lowest  resonance  of  the  galvanometer-plus-mirror  combination  is  at 
1.7  kHz.  The  convergence  time  shown  in  Fig.  5 is  8 to  10  msec,  the  design 
goal. 

The  tracking  performance  of  the  controls  is  shown  in  Fig.  6. 

An  interesting  behavior  is  observed:  the  controls  do  not  automatically  steei 

the  beam  boresight  axis  onto  the  glint  when  the  phase  and  BAT  controls  are 
used  together.  The  stable  lock  position  depends  on  the  initial  conditions 
as  shown  in  the  figures.  This  behavior  is  not  entirely  reproducible  and  was 
not  expected,  but  we  now  understand  why  it  occurs;  it  is  a consequence  of  the 
sidelobe  structure  of  the  segmented  aperture  and  of  the  choice  of  dither 
f reque  ncie  s. 

To  understand  the  multiple -state  convergence  of  the  BA  I tracking 
controls,  consider  the  diagrams  in  Fig.  7,  showing  the  element  and  array 
patterns  at  the  target  plane.  With  no  pointing/tracking  control,  the  element 
pattern  is  fixed,  and  the  COAT  sy  stem  adjusts  the  array  pattern  to  maximize 
energy  on  the  glint,  the  dotted  line  in  Fig.  7 shows  the  array  pattern  and  a 
sidelobe  for  a glint  at  point  A.  What  the  observations  in  Fig.  6 indicate  is 
that  if  a glint  is  placed  at  point  A witli  the  tracking  controls  (BAT)  off,  and 
then  the  BAT  sy  stem  is  turned  on,  the  boresight  axis  is  not  steered  to 
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Fig.  6. 

Tracking  performance  of  COAT  controls,  (a)  Tracking  control 
off.  The  formed  array  pattern  is  also  shown  for  reference. 

(b)  Tracking  control  performance.  Curves  with  tracking  control 
only  and  with  both  phase  and  tracking  controls  are  shown.  The 
ability  of  the  system  to  lock  in  a nonoptimum,  off-boresight  con- 
dition is  evident. 
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ture,  "piston-control,"  COAT  system. 
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point  A as  we  would  like.  The  array  pattern  structure  is  maintained  and  the 
BAT  system  steers  the  element  pattern  as  the  glint  moves.  This  type  of 
behavior  will  not  occur  with  a deformable  mirror  COAT  system  (or  if  the 
segmented  aperture  has  element  tilt  as  well  as  "piston"  control),  since  the 
deformable  mirror  itself  can  tilt  the  "element"  pattern  and  since  a deform- 
able mirror  does  not  produce  a side-lobe  structure  like  a segmented  array 
does. 

This  behavior  is  also  a consequence  of  the  choice  of  BAT  dither 
frequency  and  is  related  to  the  servo  loop  gain.  Intuitively,  we  can  under- 
stand what  is  happening  by  noting  that  the  servo  loop  gain  is  proportional  to 
the  change  in  target  intensity  with  respect  to  beam  motion,  and  therefore  to 
the  slope  of  the  intensity  pattern.  If  the  BAT  tracking  dither  is  fast  enough 
that  the  phase  control  cannot  follow  it,  then  the  array  pattern  is  fixed  and  it 
is  the  slope  of  the  array  intensity  pattern  that  determines  the  maximum  loop 
gain.  On  the  other  hand,  if  the  BAT  controls  dither  slowly,  the  COAT  phase 
controls  will  reform  the  array  pattern  so  that  it  is  the  slope  of  the  element 
intensity  pattern  that  is  sensed.  The  latter  case  is  what  is  necessary  for 
boresighting  to  occur. 

It  is  easy  to  show  this  result  mathematically  by  approximating  the 
total  intensity  pattern  by 


I(x)  = 


- e 


=*)' 


(i: 


where  D and  D are  the  element  and  array  peak-null  distances,  w = 1/T 
e a ’ r c c 

for  the  10%  to  90%  convergence  time  of  the  COAT  phase  controls,  and 
sine  u = sin  u/u.  By  definition,  the  COAT  system  acts  to  produce 

$ = Zttx/D  . The  BAT  tracking  dither  around  x = 0 is  x = Asinwt  and 
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sine  u = sin  u/u  = 1 - u /6. 

The  maximum  loop  gain,  G , is  found  from 
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There  are  two  cases  to  consider.  First,  when  the  BAT  dither  at  w is  faster 
than  the  phase  servo  response  (u>c),  Eqs.  (1)  and  (2)  give 


When  w<  < u>  , 
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The  ratio  of  these  two  cases  is  D^/D^,  which  is  25  for  the  DARPA/RADC 
18-element  COAT  array.  The  first  case  (Eq.  (3))  is  appropriate  for  the 
1 kHz  BAT  dither  now  in  use.  This  means  that  the  loop  gain  for  sensing  the 
element  pattern  slope  (Eq.  (4))  is  25  times  lower  than  that  which  senses  the 
array  pattern  slope  and  which  sets  the  BAT  servo  maximum  loop  gain.  The 
maximum  servo  loop  gain  for  low  frequencies  is  about  35  dB;  a factor  of 
25  lower  corresponds  to  only  7 dB.  This  small  loop  gain  is  not  enough  to 
consistently  drive  the  system  to  a boresighted  condition,  although  bore- 
sighting  will  occur  occasionally  as  indicated  in  Fig.  6. 

Use  of  a 200  to  300  Hz  dither  would  eliminate  this  behavior,  but  it 
would  also  produce  a slower  tracking  system  by  a factor  of  3 to  5.  We  have 
not  yet  demonstrated  it,  but  we  expect  that  disturbing  effects  such  as  turbu- 
lence or  defocus  will  somewhat  alleviate  the  problem  and  cause  the  BAT 
tracker  to  bo  resight  with  greater  consistency  even  with  the  1 kHz  dither. 

As  a final  result,  we  note  that  we  have  observed  glint  tracking,  black- 
hole  tracking,  and  edge  tracking  with  this  system.  Except  for  the  ineffective 
bo  re  sighting,  no  other  unexpected  behavior  has  been  observed. 


2.  Computer  Simulation  Results 

The  BAT  tracking  and  focus  controls  have  been  implemented  into  a 
multidither  computer  simulation  that  uses  a deformable  mirror  as  a dither 
and  corrector  element.  The  nonbore  sighting  effects  seen  in  the  experiment 
were  thus  not  observed.  We  also  decided  to  implement  a two-axis  focus 
control,  in  effect  using  two  cylindrical  lens  controls  that  focus  along  orthog- 
onal axes.  This  type  of  focusing  should  be  particularly  useful  with  thermal 
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blooming  since  the  negative  thermal  lens  is  much  stronger  perpendicular  to 
the  wind  direction  than  parallel  to  it.  The  computer  simulation  is  shown 
schematically  in  Fig.  8. 

The  computer  simulation  results  have  shown  the  expected  10  msec 
convergence  time  for  the  tracking  controls  and  100  msec  for  the  focus.  With 
the  focus  control,  however,  we  observed  another  initially  unexpected  result, 
the  final  converged  level  depends  on  the  magnitude  of  the  initial  defocus. 

This  behavior  is  now  understood  to  be  a generalized  manifestation  of  the 
"2mr  problem"  encountered  with  deformable  mirror  COAT  systems. 

The  behavior  of  the  focus  control  can  be  understood  by  referring  to 
Fig.  9.  Figure  9(a)  illustrates  the  well-known  behavior  of  the  on-axis 
intensity  of  a focused  beam:  as  the  focal  point  is  appic  ac,ed  from  either 
direction,  the  on-axis  intensity  goes  trough  several  maxima  and  minima. 

Each  maxima  away  from  the  focal  plane  corresponds  to  a different  Fresnel 
zone,  or  a 2mr  phase  difference  between  rays  originating  from  the  center  of 
the  focusing  lens  and  the  outer  edge.  If  the  initial  defocus  is  at  a point  such 
as  A in  Fig.  9(a),  the  BAT  focus  control  will  drive  the  focus  toward  the 
maximum  on-axis  intensity.  If  the  starting  point  is  at  B,  however,  the  BAT 
multidither  servo  moves  toward  the  nearest  maximum  at  C and  a nonoptimum 
target  intensity  results.  Figure  9(b)  shows  several  computer  runs  exhibiting 
this  type  of  behavior. 

The  results  shown  in  Fig.  9(b)  exhibit  an  unusual  behavior  when  the 
initial  defocus  is  very  large.  The  one  run  showing  several  oscillators,  and 
achieving  only  48%  convergence  at  14  psec  is  particularly  hard  to  understand. 
When  longer  run  times  for  this  case  are  used,  we  observe  that  the  conver- 
gence level  continues  to  increase,  but  very  slowly.  We  do  not  understand 
this  behavior  at  this  time,  but  we  expect  that  is  a computational  artifact 
rather  than  what  will  be  observed  in  an  experimental  system.  One  possible 
explanation  (which  has  not  been  verified  but  which  we  have  some  evidence 
for)  is  that  the  Zmr-error  sensing  and  removing  routine  built  into  the  deform- 
able mirror  simulation  may  be  cuasing  the  problem.  The  defocus  error  may 
be  large  enough  that  2u  variations  are  required  between  mirror  elements  in 
which  case  the  mirror  subroutine  is  doing  the  wrong  thing  by  removing  such 
variations.  This  point  is  further  discussed  in  Section  II-D-1. 
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Fin  9(b)  Computer  simulation  of  2-axis  BAT  focus  control  using  a 
37-element  deformable  mirror.  Computer  results  show- 
ing two  different  stable  convergence  points  which  depend 
on  the  initial  clefocus. 


We  have  made  some  computer  runs  to  determine  turbulence 
compensation  using  only  the  tracking  and  focus  controls;  no  blooming  com- 
pensation studies  have  been  made  because  of  the  difficulties  discussed  in 
Section  II-D.  Figure  10  shows  a summary  of  the  compensation  results 
achieved  with  the  BAT  controls  compared  to  those  produced  by  a 37-element 
multidither  COAT  control  system.  The  deformable  mirror  used  in  all  the 
simulation  runs  has  37  actuators.  The  assumed  propagation  scenario  is 
also  shown  in  the  figure.  For  reference,  C^Z~  21  x 10  ^ m^  ^ is  con- 
sidered to  be  strong  turbulence.  Also  shown  in  the  figure  is  the  theoretical 
limit  for  a 37-element,  deformable  mirror,  phase-conjugate  COAT  system 
that  has  perfect  phase- sensing , but  which  is  limited  by  the  ability  of  the 
mirror  to  provide  the  desired  phase  front  (see  Eq.  (22)  later  in  this  report 
for  a further  discussion  of  this  limit). 

There  are  several  noteworthy  features  of  the  data  in  Fig.  10.  First, 
the  agreement  between  analysis  and  computer  simulation  for  the  37 -element 

deformable  mirror  COAT  system  is  very  good.  The  slight  discrepancy  at 
2 

large  values  of  C^Z  may  be  caused  by  improper  operati  n of  the  simulation 
routine  designed  to  remove  2mr  ambiguities.  Second,  significant  correction 
is  obtained  with  the  tracking  or  focus  controls  when  used  separately,  but  the 
2-axis  focus  control  provides  the  most  correction.  In  fact,  when  the  track- 
ing controls  are  combined  with  the  2-axis  focus  control,  little  increase  in 

strehl  ratio  is  produced  over  that  observed  with  the  focus  control  alone.  The 

**  2 

significantly  worse  performance  of  a spherical  focus  control  above  C„Z  = 

- 1 1 1 / 3 

40  x 10  m is  also  interesting.  Finally,  not  unexpectedly,  a 37-channel 
multidither  COAT  system  produces  significantly  more  correction  than  the 
BAT  controls  alone;  the  BAT  controls  produce  roughly  a factor  of  2 improve- 
ment in  strehl  ratio  and  the  37 -channel  COAT  system  can  improve  the  strehl 
ratio  another  factor  of  2.  Experimental  data  that  are  in  good  agreement  with 
these  simulation  results  are  presented  in  Section  II-D-4. 


The  tendency  for  some  elements  of  a multidither  system  to  converge 
within  a multiple  of  2tt  of  the  desired  correction  phase  front  when  the  ini- 
tial errors  are  large. 

2-axis  focus  control  with  both  axes  driven  by  the  same  dither  and  correc- 
tion signals. 


26 


AVERAGE  GLINT  I RRADIANCE,  % OF  FREE  SPACE 


EXP  [-(0  051  h,C„,ZDT*'*)N0-,/ 


KEY  \ \ 

• NO  CORRECTION  \ \ \ 

A 2-AXIS  TILT  'V  ' 

▼ SPHERICAL  FOCUS  A.  \ \ 

□ 2-AXIS  FOCUS  \ \ 

0 2-AXIS  FOCUS  AND  TILT 

♦ 37-ELEMENT  MULTIOITHER  CONTROL 


A * 10  6Mm 
0T  « 12m 


Fig.  10.  Turbulence  compensation  produced  by  computer  simulation 
using  37-element  deformable  mirror.  Multidither  COAT 
control  is  compared  to  BAT  (tracking  and  focus  control), 
focus  control  only  (both  spherical  and  2-axis),  and  track- 
ing control  only. 


B.  COA  T/ Target-Signature  Interactions 

Contract  Amendment  No.  2 redirected  some  of  the  contract  work  to 
provide  time  and  funds  for  some  initial  studies  of  the  interaction  of  a multi- 
dither COAT  servo  system  with  backscatter  from  an  extended,  semi-diffuse 
target.  Some  of  the  results  of  these  studies  were  presented  at  the  1975 
Conference  on  Laser  Engineering  and  Applications.  ^ The  interest  in 
pursuing  these  studies  was  initially  motivated  by  the  work  of  Ogrodnik  and 
Gurski.  ^ 

1.  Problem  Statement 

The  physical  problem  known  as  "speckle  noise"  can  be  stated  briefly. 
When  a coherent  laser  beam  illuminates  a target,  some  of  the  energy  is 
scattered  back  in  the  direction  of  the  transmitted  beam.  The  backscattered 
radiation  produces  a random  intensity  pattern  (a  "speckle  pattern")  at  a 
receiver,  which  is  located  close  to  the  transmitter  for  cases  of  interest  to 
us.  When  the  appropriate  conditions  exist,  any  movement  of  this  speckle 
pattern  relative  to  the  receiver  will  produce  an  amplitude  modulation  of  the 
received  signal.  Since  a multidither  COAT  system  receiver  senses  amplitude 
modulations  on  the  backscattered  beam  that  are  produced  by  dithering  the 
transmitted  beam  phase,  spurious  amplitude  modulations  in  the  receiver  may- 
swamp  the  desired  modulations  and  thus  interfere  with  the  system  operation. 

Doppler  shifts  produced  by  rapidly  rotating  targets  may  also  produce 
false  signals,  but  these  spurious  signals  are  expected  to  be  less  important 
than  amplitude  modulation  effects  for  the  target  scenarios  of  interest.  This 
is  true  because  the  most  likely  glint  point  on  a rapidly  rotating  target  is  one 
normal  to  the  beam  incidence  directioh;  this  point  will  have  nearly  zero 
transverse  doppler  shift.  Longitudinal  doppler  shifts  are  very  narrow  band 
and  thus  will  have  minimal  effect  on  a multiple -channel  COAT  system. 

Thus  for  COAT  applications,  it  is  important  to  determine  three  tilings 
about  ta rget -signatu re  effects.  First,  what  kind  of  modulations  interfi  re  with 
the  proper  operation  of  a multidither  COAT  system?  The  interesting 
parameters  of  the  modulation  are  amplitude,  frequency  spectrum,  and  power 
spectral  density.  Second,  what  kind  of  receiver  modulations  are  produced  by 


backscatter  from  real  targets  in  expected  operational  scenarios?  That  is, 
what  scenarios  will  produce  the  modulations  that  can  affect  the  COAT 
system?  Third,  can  a multidither  COAT  system  be  designed  that  is 
insensitive  to  speckle -induced  receiver  modulation  and  if  so,  what  are  the 
critical  design  parameters?  An  expected  follow-on  contract  to  the  current 
contract  is  designed  to  produce  definitive  answers  to  all  these  questions. 
Some  preliminary  studies  of  the  magnitude  of  the  effects  are  being  pursued 
on  the  present  contract,  however. 


2. 


Experimental  Results 


For  a first  experiment,  we  chose  the  simple  extended  moving  target 
shown  in  Fig.  11(a).  The  target  is  a cylindrical  scotchlite  strip  that  could 
be  rotated  at  different  rates.  The  target  is  larger  than  the  COAT-formed 
array  pattern  in  both  dimensions  (parallel  to  and  perpendicular  to  the  rotation 
axis),  but  exceeds  the  element  pattern  diameter  only  in  one  dimension.  The 
curvature  of  the  scotchlite  provides  a localized  highlight  on  which  the  COAT 
system  can  lock.  The  laboratory  experimental  arrangement  used  is  shown 
m Fig.  11(b).  No  turbulence  or  blooming  was  introduced  in  the  paths  and 

the  BAT  controls  were  not  used.  The  receiver/transmitter  aperture  ratio 
was  chosen  to  be  0.66. 

The  experiment  consisted  of  observing  the  peak  target  irradiance  and 
the  COAT  receiver  signal  frequency  spectrum  when  the  target  was  stationary 
and  moving.  Figure  12  shows  the  observed  receiver  power  spectra  for 
several  cases.  The  figures  show  a spectrum  analyzer  output  when  the  input 
is  the  COAT  photomultiplier  signal.  The  spectrum  analyzer  bandwidth  was 
10  Hz  in  Figs.  12(a)  and  (b),  and  200  Hz  in  (c)  and  (d).  Each  photograph  con- 
sists of  five  2-minute  scans  over  a 0 to  50  kHz  range.  Figures  12(a)  and  (c) 
are  the  observed  spectrum  when  the  target  is  stationary  and  Figs.  12(b)  and 
(d)  are  the  spectrum  when  the  target  rotates.  The  target  rotation  speed  was 
chosen  to  give  the  largest  observed  change  in  the  spectrum  over  the  dither 
band.  For  the  stationary  targets,  the  noise  level  is  low  and  the  dither  signals 
are  clearly  evident  above  the  noise.  When  the  target  is  moving,  the  noise 
level  increases  significantly  and  most  of  the  dither  signals  are  no  longer 
evident.  The  peak  target  irradiance  was  the  same,  however,  whether  the 
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Fig.  11.  Experimental  apparatus  for  COAT/target  studies,  (a)  Rotating 
scotchlite  cylinder  used  for  target,  (b)  Schematic  of  optical 
setup. 
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Fig.  12. 


target  was  stationary  or  moving.  For  this  simple  target,  if  the  COAT 
system  could  converge  on  the  target  when  it  was  stationary,  the  system 
could  converge  equally  well  when  the  target  was  moving.  For  every  target 
and  target  motion  rate  that  we  have  tried  to  date,  this  same  observation  has 
been  made:  no  noticeable  effects  when  the  target  moves. 

The  observation  of  no  effect  on  the  COAT  system  performance  is 
contrary  to  a statement  made  in  the  third  management  report.  The  very 
earliest  experiments  did  seem  to  show  some  effect.  What  we  found,  how- 
ever, was  that  very  large  extraneous  dither  signals  were  being  introduced 
into  the  COAT  receiver  by  ground  loops,  not  by  target  speckle  effects.  Since 
the  ground  loops  have  been  eliminated,  we  have  not  been  able  to  find  a target 
which  "confused"  the  COAT  system  solely  because  of  its  motion. 

Since  we  have  seen  no  deleterious  COAT-system/target-modulation 
effects,  a natural  question  is:  "What  kind  of  spurious  signals  can  reduce  the 

system  performance?"  We  thus  undertook  an  experiment  whose  goal  was  to 
generate  large  amplitude  spurious  amplitude -modulated  signals  capable  of 
overloading  the  dither  servo  channels.  In  this  experiment,  the  target  was  a 
single  glint  and  a large  receiver  aperture  was  used.  Spurious  amplitude 
modulations  were  artificially  introduced  into  the  COAT  system  by  putting  an 
acousto-optic  modulator  (AOM)  in  the  optical  beam  path,  as  shown  in  Fig.  13. 
The  AOM  amplitude  modulates  the  transmitted  beam  at  a frequency  set  by  ' 
the  voltage -controlled-oscillator  (VCO),  shown  in  F’g.  13.  Figure  14  shows 
the  modulated  light  incident  on  the  target  glint  when  the  oscillator  was  not 
swept  in  frequency.  Modulation  depths  of  64%  were  observed  over  the  dither 
frequency  band  (8  to  32  kHz)  with  a triangular -wave  drive  to  the  VCO. 

Manually  tuning  the  VCO  through  tire  dither  passband  sequentially  over- 
loaded the  feedback  loop  at  each  dither  frequenc>  and  caused  a 10%  reduction 
in  peak  target  ir radiance  (a  loss  of  proper  phasing  in  1 of  the  18  COAT 
channels).  Spectrum  analyzer  measurements  indicated  that  the  spurious 
signal  level  of  100  mV  generated  in  this  manner  was  much  larger  than  the 
3 mV  dither  amplitudes.  This  initial  result  establishes  the  fact  that  external 
modulation  applied  to  the  COAT  beam  could  overload  the  control  circuitry 
in  any  single  channel  and  cause  the  loss  of  that  channel  in  forming  the 
converged  beam. 
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Fig.  13.  Experimental  setup  to  determine  effect  of  spurious  amplitude 
modulations  in  a COAT  receiver. 
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Fig.  14.  Light  intensity  incident  on  target  with  acousto- 
optic modulator  (AOM)  in  the  optical  path. 


To  roughly  simulate  the  effect  of  worst-case  broadband  speckle 
modulations,  the  VCO  was  swept  over  a selectable  bandwidth  with  sweep 
periods  of  1 msec  and  0.  5 msec  (shorter  time  than  the  COAT  system  con- 
vergence time  of  1.  5 to  2.  0 msec).  During  this  series  of  tests,  the  COAT 
receiver  frequency  spectrum  was  recorded  and  the  power  level  in  the  con- 
verged beam  was  measured  with  the  sweep  voltage  turned  on  and  off.  The 
resulting  spectra  for  six  different  sweep  bandwidths  and  a 1 msec  sweep 
period  (0.5  msec  across  the  sweep  bandwidth)  are  shown  in  Fig.  15.  The 
center  frequency,  f , was  selected  for  maximum  effect  on  the  COAT  system. 
The  sweep  bandwidths  in  Fig.  15  range  from  1 kHz  to  30  kHz  (f  = f ± BW/2), 
but  the  ratio  of  the  peak  target  ir  radiance  with  the  spurious  modulation  to 

that  without  it  (P  ,/P  ) is  nearly  constant.  A value  of  P ,/P  =0.69  = 

^ mod  o 7 mod  o 

(15/18)  corresponds  to  a loss  of  proper  phasing  in  three  channels  on  the 
average  and  (14/18)  = 0.60  corresponds  to  a loss  of  4 channels.  Similar 

results  are  obtained  when  the  sweep  period  is  reduced  to  0.  5 nsec.  A graph 
of  the  number  of  active  elements  lost  versus  modulation  bandwidth  is  shown 
in  Fig.  16(a)  and  the  corresponding  loss  of  converged  beam  power  loss  versus 
modulation  bandwidth  is  shown  in  Fig.  16(b). 

The  power  spectrum  of  the  triangular-wave  swept  modulation  used 
here  consists  of  a sum  of  all  the  harmonics  of  the  sweep  frequency  spaced 
about  the  center  frequency: 


00 

P = ^ a exp[*i27r  (f  - nf  )"|  (5 

m n L ° sweep/J 

n = - oo 

where  the  a^  are  determined  in  part  by  the  sweep  bandwidth.  Although  this 

type  of  power  spectrum  will  probably  not  be  characteristic  of  any  target- 

induced  speckle  spectrum,  it  does  have  several  properties  which  a speckle- 

induced  receiver  spectrum  should  also  have.  First,  the  bandwidth  of  the 

extraneous  receiver  signals  is  limited,  but  will  almost  certainly  exceed 

4 

35  kHz  in  cases  of  most  interest.  Second,  the  product  of  peak  power 
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MODULATION  BANDWIDTH,  kHz 


(a)  Observed  peak  target  irradiance. 

Effect  of  acousto  optic  modulation  on  DARPA/RADC 
COAT  system  as  a function  of  total  modulation 
bandwidth. 


spectral  amplitude  and  bandwidth  is  constant.  The  experimental  AOM  data 
obey  such  a relationship  as  demonstrated  in  Fig.  17.  What  this  means  is 
that  the  receiver  can  see  a 100%  average  modulation-depth  signal  at  only  a 
single  frequency;  as  the  modulation  bandwidth  increases,  the  average  receiver 
voltage  modulation  depth  (not  the  instantaneous  depth  or  amplitude)  goes  down 
as  th  reciprocal  of  the  square  root  of  the  bandwidth. 

Determination  of  the  actual  spectrum  of  a real  speckle -induced 
receiver  signal  is  beyond  the  scope  of  this  contract.  As  a worst  case,  how- 
ever, consider  a spectrum  of  equal  amplitude  sinusoidal  frequencies,  spaced 
by  100  Hz  over  a 35  kHz  bandwidth,  which  is  close  to  a worst-case  situation. 
Since  there  are  350  frequencies,  the  average  receiver  voltage  modulation 
depth  for  any  one  must  be  (350)“^  ^ = 0.  042,  or  4.  2%  to  maintain  a constant 
power-bandwidth  product.  This  assumes  a 100%  speckle  modulation  depth, 
which  will  occur  only  for  purely  diffuse  targets.  If  we  take  a 50%  modulation 
depth  since  realistic  targets  will  probably  have  some  specular  returns,  the 
average  modulation  depth  drops  to  2.  1%.  This  value  is  below  the  5.6%  peak 
receiv  ; r modulation  produced  by  a multidither  system  operating  with  a 
±20°  dither  modulation. 

A further  consideration  that  we  have  not  taken  into  account  is  that  the 
spectral  components  of  speckle -induced  receiver  noise  will  have  randomly 
varying  phases  with  respect  to  the  dither  modulation.  The  servo  synchronous 
detectors  will  thus  further  discriminate  against  this  type  of  noise.  We 
expect  to  have  an  analytical  demonstration  of  the  magnitude  of  this  discrimi- 
nation sometime  in  the  near  future. 

These  brief  studies  seem  to  indicate  that  target -signature  effects  are 
not  going  to  be  as  severe  as  was  thought  earlier.  More  realistic  targets 
and  motions  need  to  be  studied,  however,  and  other  effects  such  as  more 
servo  channels  or  receiver  signal-to-shot  noise  need  to  be  investigated.  The 
presence  of  turbulence  or  blooming  may  also  have  some  effect. 


A more  realistic  speckle -induced  spectrum  might  be  a uniform  triangular 
or  Gaussian  distribution  with  a 50  kHz  bandwidth. 
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3. 


Analytical  Results 


During  the  last  quarter  the  multidither  servo  computer  simulation 
developed  at  HRL  was  modified  to  include  the  effects  of  speckle  modulation 
on  system  convergence.  The  results  of  the  simulation  have  also  been 
compared  to  a statistical  analysis  that  treats  both  the  COAT  system  and  a 
model  of  the  speckle  noise. 

The  computer  simulation  was  modified,  as  shown  in  Fig.  18,  to 
include  a multiplic~five  receiver  noise  — the  speckle -induced  noise.  The 
amplitude  modulation  noise  chosen  for  these  initial  studies  has  the  following 
form: 


M 

s 


N M 

7 + i 1 ajk  sin  ( V + V 
k=l  j=l 


(6) 


whe  re: 


N 

M 


wjk 

V 

ajk 


number  of  channels 

number  of  erroneous  signals  near  each  dither  frequency 

the  jtii  frequency  in  a band  of  M equally  incremented 
frequencies  centered  at  the  dither  frequency 

a randomly  chosen  phase  angle  between  0 and  2n 

a randomly  chosen  amplitude  between  0 and  an  adjustable 
maximum,  a 

max 


A topical  spectral  breakdown  of  this  noise  is  shown  in  the  inset  to  Fig.  19(c). 
There  are  21  discrete,  equally -spaced  frequencies,  at  or  near  each  dither 
frequency  of  the  COAT  system,  whose  maximum  amplitude  is  set  by  the 
parameter  a . Each  noise  frequency  component  has  a randomly  chosen 
amplitude  and  phase.  This  signal  multiplies  the  usual  COAT  receiver  signal. 
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Fig.  19.  COAT  computer  simulation  time  histories. 

(a)  Broadband  shot  noise  signal-to-noise 
ratio  = (S/N shot  = 1010  and  a = 0 (no 
speckle  noise),  (b)  a = 0.  008  corres- 
ponding to  <MS>  = 0.  064.  (c)  Speckle 

noise  modulation  used  to  produce  (b). 


For  the  simulation  studies  performed  to  date,  the  21  noise  frequencies 
are  spaced  ±8  3 Hz  around  each  dither  frequency.  This  is  a rather  arbitrary 
choice,  but  it  does  insure  that  all  spurious  modulation  is  well  within  the 
400  Hz  servo  bandwidth.  An  amplitude  limitation  was  also  included  to  keep 
the  return  signal  from  exceeding  an  absolute  value  of  unity  just  before 
entering  the  loop  gain  amplifier  (noted  as  Gj  in  Fig.  18).  This  modification, 
which  avoids  instabilities  resulting  from  the  automatic  gain  control  mech- 
anization, is  used  in  the  lab  model  of  the  COAT  array.  Simulations  were  run 
at  three  values  of  broad-band  signal -to -noise  ratio  (receiver  shot  noise): 

1010,  40,  and  20.  When  (S/N)  h = 1010  and  amax  = 0 (no  speckle  noise), 
the  simulation  produces  a time  history  of  the  power  on  the  target  glint  as 
shown  in  Fig.  19(a).  When  the  speckle  noise  modulation  shown  in  Fig.  19(c) 
is  present  (a  = 0.008  and  <M  > = 0.064  - See  Eq.  (6)),  the  time  history 

is  pieseiu.  Vamax  s 

in  Fig.  19(b)  results.  Notice  that  although  the  rms  value  of  Mg  is  only  0.  064 

(around  M =0.  50),  the  maximum  peak-to-peak  value  is  about  0.4. 

' s 

The  behavior  shown  in  Fig.  19(b)  for  amax  ^ 0 is  typical.  In  general, 
the  speckle  noise  has  two  effects:  (1)  the  maximum  percent  convergence  level 

is  lowered;  and  (2)  a random  oscillation  of  the  intensity  occurs  between  the 
new  maximum  convergence  level  and  a minimum  level.  This  oscillation  is 
much  slower  than  the  dither  oscillations  and  its  magnitude  varies  with  con- 
vergence level.  Figure  20  shows  a plot  of  the  maxima  and  minima  recorded 
during  these  simulations  as  a function  of  amax-  The  °Pen  looP  &ain>  Gl>  was 
the  same  for  each  case.  At  high  average  convergence  levels  (low  values  of 

a ),  the  random  oscillations  have  much  smaller  magnitudes  than  at  the 
max 

low  convergence  levels  (high  values  of  amax)- 

In  order  to  understand  the  behavior  of  the  multidither  COAT  system 

with  multiplicative  speckle  noise,  we  have  begun  a statistical  analysis  of 
such  a COAT  system  which  includes  the  speckle  noise.  The  analysis  is  not 
yet  complete,  but  it  does  provide  some  understanding  of  the  COAT  system 
ope  ration. 

The  object  of  the  analysis  is  to  explain,  and  if  possible,  to  predict 
the  behavior  observed  in  the  computer  simulations.  Note  that  the  voltage 
generated  at  the  photomultiplier  is  proportional  to  the  intensity  of  the  radia- 
tion incident  on  it.  The  servo  error  information  contained  in  this  signal 
consists  of  modulations  at  the  dither  frequencies,  u>n.  If  the  reflected 
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intensity  is  further  modulated  by  the  speckle  noise  function  Mg,  false  error 
information  is  generated  in  the  servo  loop.  The  voltage  at  the  dither  f-e- 
quencies  will  then  consist  primarily  of  the  sum  of  the  correct  error  infor- 
mation plus  the  false  error  signal.  Other  sources  of  noise  at  these  fre- 
quencies are  ignored  in  this  analysis. 

As  we  have  just  stated,  the  total  error  signal,  S , is  just  the  sum  of 
Sjy  the  correct  dither  error  signal,  and  Sg,  the  false  error  signal  due  to 
speckle  (both  in  volts).  The  strength  of  the  dither  signal  depends  on  the 
convergence  level  while  the  strength  of  the  speckle  signal  depends  on  both 
the  convergence  level  and  the  parameter  a . Based  on  preliminary  cal- 
culations,  we  hypothesize  that  the  time -averaged  convergence  level  main- 
tained by  a COAT  system  with  a speckle  modulation  M^,  corresponds  to  a 
fixed  ratio  of  the  average  dither  signal  power  to  the  average  speckle  modula- 
tion signal  power.  With  this  hypothesis  in  mind,  we  will  calculate  the 
expectation  values  of  these  average  powers  as  a function  of  the  parameters 

a and  the  mean  convergence  level  I.,,  which  will  be  defined  shortly, 
max  M 1 

Expressions  will  be  derived  which  will  enable  us  to  relate  the  three  param- 

2 2 

eters  a , I.,,  and  the  power  ratio  p = S^/S-.  . This  will  allow  us  to  cal- 
max  M D S 

culate  I as  a function  of  a at  fixed  values  of  p.  If  the  above  hypothesis 
M max  r ’ r 

is  correct,  then  these  curves  should  predict  the  same  kind  of  behavior  as 
the  simulation  data  in  Fig.  ZO. 

Consider  the  electric  field  of  N separate  waves  having  the  same 
amplitude  and  wavelength  incident  on  a given  point.  Assume  the  polariza- 
tion of  each  wave  is  adjusted  so  that  they  are  all  plane  polarized  parallel  to 
each  other.  Then  the  total  field,  E^,  is 


irJt) 


ET(t)  = LEoC 


where  Eq  is  the  magnitude  of  each  wave  (assumed  equal  for  each  wave),  and 
T (t)  is  the  timc-\arying  phase  of  the  n*^  wave  at  the  target  glint. 

The  radiant  intensity  is  equal  to  the  absolute  square  of  the  electric 
field.  Consequently,  a radiant  intensity  normalized  to  unity  when  all  T - 0 
may  be  written  as: 


(8) 
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In  the  COAT  multidither  array  the  N phase  angles  are  physically^ 
dithered  at  different  frequencies.  For  such  an  array  the  phase  in  the  n 
channel  T (t)  = (3  (t)  + + sin  c^t.  In  this  situation  ^(t)  represents  a mean 
phase  ang’ie  and  4 is  the  amplitude  of  the  small  dither  deviation  from  the 
mean  (typically  * < tt/6  radians).  Substituting  TJt)  = P„(t)  + - sin  c^t  into 
eq.  (8)  and  expanding  the  result  in  Bessel  functions  gives 
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Since  we  are  only  interested  in  the  dither  error  signal  and  the  convergence 
level,  we  need  to  consider  only  the  first  Uvo  terms  in  the  Bessel 

expansions: 
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Thus  we  car.  separate  w... 
modulations  at  the  dither  frequencies 
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te  out  I , the  mean  convergence  level,  and  1^,  the 
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We  can  now  express  the  intensity  on  the  glint  as 


ZT  (IM  + ZD  + V ZM  + *D  ’ 


(12) 


where  1^  is  the  sum  of  all  modulations  in  eq.  (9)  that  are  not  included  in 
Ip  (i.  e.  , modulations  at  harmonics  plus  sum  and  differences  of  the  dither 
frequencies).  These  can  be  considered  as  noise  since  they  are  not  used  in 
the  servo  system  to  generate  error  information  and  since  they  are  higher 
order  terms  in  the  variable  *\>. 

In  the  simulations,  the  total  intensity  incident  on  the  photomultiplier 
is  given  by  the  product  M^I^,.  Combining  eqs.  (6)  and  (12)  gives 
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The  second  and  third  terms  in  eq.  (13)  are  the  ones  we  are  interested 
in.  The  second  term  contains  the  correct  dither  information,  while  the 
third  contains  the  false  dither  signal  from  the  speckle  noise.  The  first 
term  is  a dc  term  representing  the  mean  intensity  at  the  receiver  and  thus 
is  not  important  to  our  discussions.  The  last  term  consists  of  intensity 
modulations  at  or  near  harmonics  and  sums  and  differences  of  dither 


frequencies.  This  term  can  be  neglected  for  two  reasons.  First,  most  of 

the  modulations  will  be  outside  the  dither  band  and  will  therefore  be  filtered 

out  by  the  low-pass  filter  in  the  servo  electronics.  Second,  since  it  is 

proportional  to  1^,  it  is  an  order  of  magnitude  smaller  than  the  third  term, 

which  is  proportional  to  I . 

M 

The  voltages  and  Sg  are  thus  proportional  to  the  second  and 
third  terms  respectively.  Letting  K be  the  photomultiplier  proportionality 
constant,  we  have 
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When  the  COAT  system  reaches  a steady  state  convergence  level 

(such  as  that  shown  in  Fig.  19(b)  for  example),  we  can  calculate  the  average 

dither  signal  power,  P , and  the  average  speckle  signal  power,  P , 

LJ  5 

normalized  to  unit  resistance  by  using  eqs.  (14)  and  (15).  The  results  are 
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We  now  want  to  calculate  the  power  ratio  as  a function  of  convergence 
level  or,  in  other  words,  as  a function  of  I In  order  to  do  this,  we  need 
a statistical  model  of  an  N-channel  COAT  array  at  a steady  state  conver- 
gence level.  This  will  allow  us  to  calculate  expectation  values  of  P . P 
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and  any  other  desirable  quantities  as  a function  of  a statistical  parameter 
a,  which  characterizes  the  convergence  level. 

In  order  to  define  the  parameter  a,  we  first  recognize  that  the  phase 
angles  (3^  are  distributed  over  a limited  overall  angular  spread.  We  define 
the  magnitude  of  this  spread  as  a.  This  can  best  be  understood  in  terms  of 
a phasor  diagram.  In  Fig.  21,  for  example,  we  show  the  N phases  as  more 
or  less  uniformly  distributed  over  an  angle  a.  As  a approaches  zero,  the 
laser  array  will  converge  (i.  e.  , 1^,  -♦  1 and  thus  I 1).  Conversely,  as 
a approaches  2ir,  the  array  approaches  a convergence  level  of  zero.  The 
range  of  all  possible  relative  phase  differences,  (3^  - (3  , between  any  two 
channels,  at  a specified  convergence  level,  is  thus  0 to  a.  The  quantity  a 
therefore  characterizes  the  convergence  level,  recognizing,  of  course,  that 
this  is  true  only  in  a statistical  sense.  We  now  make  the  rather  arbitrary 
assumption  that  for  a given  convergence  level,  the  N phases  (3  are  uni- 
formly distributed  between  a and  a + a as  shown  in  Fig.  21.  We  do  not 
expect  our  results  to  depend  strongly  on  the  assumed  distribution,  and  a uni- 
form distribution  simplifies  the  computations.  The  probability  density  that 
a given  (3^  = (3  is  just  P((3)  = a * over  this  range  and  zero  everywhere  else. 
The  probabiliy  density  P(Ap),  for  an  arbitrary  (3^  - (3^  = A(3,  is  the  auto- 
correlation function  of  P((3),  given  by 


P(A(3) 


= i . 

a 2 

Q 


(18) 


This  expression  holds  in  the  region  -a  < Ap<  a.  Outside  that  region  P(Ap) 

h 0. 

Knowing  the  probability  distribution  allows  us  to  calculate  various 
expectation  values.  Consider,  for  example,  the  expectation  value  of  I . 
Using  the  definition  of  I in  eq.  (10)  gives 


1 J>>  N N 
<im>  -i ? + -rr 


N 


X X <CO$  ,|3n  - 


(19) 


n- 1 1 

n/m 


Using  eq.  (18),  it  is  easy  to  show  that 
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(20) 


<cos(0n  - Pm)>  = f P(A(3)  cos(Ap)d(A|3)_  = ( 1 - cosa)  . 

J a 

-co 


Thus  we  have 


<IM> 


b + Jo,w  1 1 


('  - i)  5 


( 1 - cos  a) 


(21) 


Figure  22  shows  a plot  of  versus  a for  an  18-element  array,  and  a 

peak  dither  amplitude  of  4*  = 0.  349  radians  (20°). 

During  the  next  quarter  we  will  use  this  analysis  to  calculate  expec- 
tation values  of  P^,  P and  their  ratio,  p,  as  functions  of  I.,.  Plots  of  I 

US  MM 

versus  a at  constant  power  ratio  will  then  be  made  to  compare  with 
max  r r 

simulation  data  such  as  that  in  Fig.  20. 
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Fig.  21.  Phaser  diagram  of  a partially 

converged,  N-element  COAT 

array  defining  the  quantities  a 

and  a used  in  the  text, 
o 


c. 


Thermal  Blooming  and  Turbulence  Compensation  Experiments 

1.  Artificial  Turbulence  Generation 

As  discussed  in  the  previous  report,  * we  have  chosen  to  use  a single 
phase  screen  to  produce  turbule nee -like  effects.  We  have  produced  tv/o 
phase  screens  by  sputtering  SiO^  onto  a thin  quartz  plate.  A scan  across 
one  of  these  plates  made  by  a Dektak  electronic  micrometer  is  shown  in 
Fig.  23.  The  peak-to-peak  optical  phase  variation  is  about  3/4  wavelength 
and  the  minimum  size  of  a given  "turbule"  is  about  0.  5 mm.  For  significant 
turbulence  effects,  this  diameter  should  be  greater  than  the  width  of  a single 
element  in  the  COAT  transmitted  array,  but  less  than  the  full  beam  diameter 
(five  times  one  element  width). 

Because  of  the  optical  arrangement  that  we  have  been  using,  these 
plates  have  not  produced  enough  ‘'um  distortion.  The  mounting  technique 
produces  significant  beam  steering,  but  beam  breakup  and  scintillation  is 
minimal.  During  the  next  quarter  we  will  fabricate  plates  with  "stronger" 
turbulence  and  improve  the  optical  arrangement  to  maximize  the  effect. 

For  the  experiments  discussed  below,  we  have  placed  saran  wrap  over  one 
of  the  quartz  plates.  This  gives  a more  severe  distortion  — sometimes  far 
too  much  — but  also  produces  a lot  of  scattering  or  amplitude  effects.  This 
arrangement  is  being  used  only  for  our  p eliminary  tests. 

2.  Flowing  Gas  Cell 

The  flowing  gas  cell  discussed  briefly  in  the  previous  report  has  been 
constructed  and  tested.  The  initial  construction  had  several  leaks  which  have 
now  been  eliminated.  The  unit,  shown  in  Fig.  24,  can  produce  turbulence  - 
free  air  flow  of  up  to  12  cm/sec  as  measured  by  a hot-wire  anemometer; 
similar  performance  is  observed  with  a mixture  of  Xenon  and  NO^. 

We  have  found  it  necessary  to  establish  the  absorption  of  this  cell 
by  direct  periodic  measurement  rather  than  by  filling  it  with  a fixed  NO^ 
partial  pressure.  This  is  necessary  because  the  absorption  changes  slowly 
and  in  an  unpredictable  manner  as  the  cell  runs.  The  cell  transmission  will 
vary  as  much  as  30%  over  a two-hour  period  in  an  oscillatory  manner.  We 
do  not  understand  this  behavior,  but  attribute  part  of  it  to  cleanup  and  then 
re-emission  of  the  NO^  by  the  sintered  nickel  diffuser  plates  in  the  cell. 
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(a)  Schematic  drawing.  Sintered 
nickel  plates  are  used  for  dif- 
fusers to  produce  laminar 
flow. 


FLEXIBLE- 

STAINLESS 

TUBING 


OPTICAL 

PATH 


BLOWER 


VALVE 


VALVE 


(b)  Cell  is  mounted  with  a vertical  path  to 
minimize  buoyancy  effects. 


Fig.  24.  Flowing  gas  cell  or  blooming  studies 


3.  Blooming  Compensation 

Because  of  the  initial  difficulties  with  the  flowing  gas  cell,  we  have 
been  able  to  make  only  a few  preliminary  measurements  of  blooming  com- 
pensation. We  have  reworked  the  optical  arrangement  now  so  that  a 62  cm 
focal  length  lens  focuses  the  beam  through  the  cell.  This  puts  the  blooming 

in  the  last  81%  of  the  propagation  path.  The  observed  results  are  shown  in 

15  6 

Fig.  25;  previous  data  ’ ’ are  shown  in  Fig.  26  for  comparison. 

The  data  in  Fig.  25  were  obtained  using  the  18  channels  of  phase 
control  and  the  tracking  controls.  The  tracking  controls  alone  produced  no 
observable  change  in  the  peak  target  irradiance.  They  could,  of  course, 
increase  the  peak  power  on  a selected  glint  point  by  steering  the  beam  onto 
that  point.  The  improvement  in  peak  target  irradiance  is  about  35%,  still 
far  below  the  factors  of  2 to  3 observed  with  predictive  phase  compensation. 

We  still  do  not  understand  why  we  see  no  significant  compensation. 

It  is  apparent  from  the  data  in  Fig.  26(a)  that  blooming  distortions  in  the 
transmitter  near  field  can  be  removed;  it  is  equally  apparent  from 
Fig.  26(b)  that  those  in  the  far  field  cannot  be  removed.  It  is  difficult  to 
believe  that  it  is  only  the  blooming  in  the  initial  19%  of  the  propagation  path 
that  can  be  compensated,  but  the  data  in  Figs.  25  and  26  ten!  to  point  that 
way.  During  the  next  quarter  we  will  look  at  two  cases  still  of  interest: 

(1)  blooming  in  the  first  81%  of  the  propagation  path  (as  opposed  to  the  last 
81%  as  in  Fig.  25);  (2)  slewing  beams.  The  second  case  is  possible  now 
with  the  tracking  controls  operational. 

4.  Turbulence  Compensation 

We  have  made  some  turbulence  compensation  studies  using  the 
artificial  phase  screen  discussed  in  Section  II-C-1  above.  We  have  looked 
at  four  cases:  (1)  NO  COAT  control;  (2)  phase  control  only;  (3)  tracking 

control  only;  (4)  phase  and  tracking  control.  The  experimental  arrangement 
is  shown  in  Fig.  11(b)  and  Table  I summarizes  the  observations.  As  can  be 
seen  from  Table  I,  the  phase  and  steering  distortions  in  our  artificial 
turbulence  have  roughly  equal  contributions  to  the  uncompensated  distortion. 
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PEAK  TARGET  IRRADIANCE, 
ARBITRARY  UNITS 


4184  11 


INPUT  POWER.  ARBITRARY  UNITS 


Fig.  25.  Thermal  blooming  compensation  data  obtained  with  flowing  gas 
cell.  Also  shown  is  the  focused  propagation  path  relative  to 
the  blooming  medium. 
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b)  FAR  FIELD  BLOOMING 


Fig.  26.  Thermal  blooming  compensation  for  two  propagation  scenarios 
different  from  that  in  Fig.  25.  (a)  Thin  static  liquid  cell  in 

near -field  of  the  transmitter,  (b)  Static  gas  cell  (but  trans- 
verse wind  produced  by  mechanical  motion)  in  last  55%  of 
propagation  path. 
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Table  I.  •'-COAT  Compensation  for  Artificial  Turbulence 


Case 

Peak  Target 
Irradiance, 
arb.  units 

Strehl  Ratio, 
S 

NO  COAT 

0.  39 

0.  65 

COAT  phase  control 

0.  52 

0.  87 

COAT  tracking  control 

0.  52 

0.  87 

Phase  and  tracking 

0.  57 

0.  95 

control 

S = 1 is  defined  for  no  turbulence,  COAT  phase  control. 

D.  Computer  Simulation  Results 

1.  Turbulence  Compensation 

9 10 

We  have  demonstrated  elsewhere  ’ that  the  effect  of  atmospheric 
turbulence  on  the  ensemble -ave raged  irradiance  distribution  can  be  accounted 
for  by  introducing  a random  phase  screen  in  front  of  the  transmitting  aperture. 
The  appropriate  phase  distribution  is  that  which  would  be  produced  by  a point 
source  located  on  the  target.  This  implies  that  an  ideal  adaptive  optics  sys- 
tem should  be  able  to  completely  compensate  for  the  effects  of  the  turbulence 
by  introducing  the  conjugate  phase  distribution  as  a correction.^  The  degree 
to  which  this  is  achieved  depends  on  whether  the  target  has  a point-glint 
scatterer  that  provides  the  requisite  phase  information  and  on  the  degree  to 
which  the  adaptive  optics  system  can  reproduce  the  desired  phase  distribution. 


In  this  discussion,  we  assume  that  amplitude  scintillation  effects  are  small. 
If  they  are  not  small,  the  effect  of  atmosphere  turbulence  cannot  be 
represented  simply  by  a phase  screen  but  rather  we  must  also  introduce 
a random  apodization. 

^Subject  to  the  assumption  that  amplitude  effects  are  negligible.  t 
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we 


Previously*^’  10  we  have  demonstrated  that  phase  conjugate  COAT 

control  with  a deformable  mirror  can  correct  a beam  focused  through  a 

turbulent  atmosphere  back  to  within  70%  to  95%  of  the  diffraction -limited 

.10 

(free-space)  target  irradiance.  We  have  also  demonstrated  analytically 

that  a multidither  COAT  system  will  arrive  at  the  same  phase  correction  as 

a phase  conjugate  system  when  correcting  for  linear  (not  power -dependent) 

distortions.  The  amount  of  correction  depends  on  the  range  (Z),  turbulence 

strength  (C^T),  wavelength  (\),  transmitter  diameter  (DT),  and  number  of 

deformable  corrector  mirror  actuators  elements  (N  ).  The  dependence  of 

• J ..  . . 10 

the  residual  phase  error  after  correction  is  described  by 


> 


Z D 


5/3 


(22) 


where  <d>  is  the  residual  mean  square  phase  error  across  the  aperture 
e 

after  correction.  This  result  includes  an  approximation  for  the  shape  of  a 
deformable  mirror  surface.  The  effective  strehl  ratio  is  given  by 


S ~ exp  <4>e2>J 

The  value  of  S is  plotted  in  Fig.  28  versus  C^  for  several  cases  of  interest. 

To  provide  a check  on  the  above  analytical  work,  we  have  performed 
two  types  of  mirror  simulations.  In  each,  the  effect  of  a turbulent  phase 
screen  placed  in  front  of  an  aperture  was  compensated  by  a mirror  with  a 
finite  number  of  actuators.  The  average  Strehl  ratio  obtained  from  a sequence 
of  five  independent  phase  screens  was  determined  as  a function  of  the  number 
of  actuators.  In  one  of  the  simulations  the  deformable  mirror  software  was 
used  in  conjunction  with  a sinusoidal  multidither  COAT  algorithm  of  the  type 
used  previously  at  HRL  in  our  COAT  servo  system  studies.  In  the  othei 
simulation,  the  mirror  was  modeled  by  a segmented  mirror  with  piston  and 
tilt  control  on  each  segment.  The  piston  and  tilt  settings  for  each  segment 
were  determined  by  a least  squares  fit  to  the  random  phase  surface  over  the 
segment. 
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= STREHL  RATIO 


Fig.  27.  Strehl  ratio  associated  with  atmospheric  turbulence 
sus  number  of  mirror  actuators. 


The  agreement  between  the  piston  and  tilt  results  and  the  theoretical 
predictions  obtained  from  Eqs.  (22)  and  (23)  is  very  good  (see  Ref.  9 and 
Fig.  10  of  this  report),  which  gives  us  confidence  in  the  theoretical  results 
given  in  these  equations.  We  believe  that  these  results  are  representative 
of  those  that  would  be  obtained  with  a deformable  mirror  in  the  absence  of 
the  2Ntt  problem  discussed  below.  The  results  actuall  obtained  with  a de- 
formable mirror  will  depend  on  the  degree  to  which  the  2Ntt  problem  is 
avoided. 

The  agreement  between  the  theory  and  the  deformable  mirror- 

multidither  COAT  simulation  results  is  reasonably  good  for  moderate  values 

2 2 

of  but  is  poor  for  large  values  of  C^.  We  attribute  the  poor  results 
obtained  at  large  values  of  to  a 2NTi-type  of  behavior.  A similar  problem 
is  observed  with  the  focus  control  computer  simulation  (see  Fig.  9(b)).  The 
deformable  mirror  simulation  that  was  used  in  these  runs  has  a 2Ntt  correc- 
tion loop  that  introduces  a 2tt  correction  whenever  the  phase  difference 
between  actuators  exceeds  four  radians.  The  intent  is  to  suppress  2ir  errors 
introduced  by  the  servo  system.  However,  if  the  phase  distortion  that  is 
to  be  removed  by  the  mirror  changes  by  more  than  four  radians  between 
actuators,  the  2t r "correction"  introduced  by  the  2Ntt  correction  loop  is,  in 
fact,  a 2tt  error  and  the  mirror  performance  is  corresponding  degraded. 

This  problem  could  be  avoided  by  removing  the  2Nit  correction  loop  but  then 
we  would  be  faced  w'ith  2Ntt  servo  errors.  One  way  to  avoid  this  problem  is 
to  use  more  actuators  so  that  the  phase  change  between  actuators  never  ex- 
ceeds four  radians.  Another  way  would  be  to  design  the  mirror  so  that  the 
likelihood  of  2Ntt  errors  is  reduced,  in  which  case  the  2Ntt  correction  loop 
could  be  eliminated.  More  work  clearly  needs  to  be  done  on  this  problem. 

2.  Blooming  Compensation 

Although  a multidither  system  and  a phase  conjugate  system  arrive 
at  the  same  answer  for  linear  propagation  distortions,  the\  do  not  reach  the 
same  result  for  nonlinear  thermal  blooming  distortions.  In  the  case  of 
blooming,  it  can  be  demonstrated  that  applying  the  phase  conjugate  ot  a wave 
generated  by  a point  source,  that  propagates  from  the  focal  plane  back  to  the 
transmitter  does  not  insure  maximum  target  irradiance: 


tin*  correction 


algorithm  is  not  equivalent  to  a true  target  irradiance  maximization 
algorithm.  A multidither  algorithm,  however,  is  a true  maximization 
scheme  (barring  possible  problems  with  2tt  ambiguit.es). 

In  trying  to  simulate  multidither  compensation  for  blooming,  we  have 
been  stopped  by  the  lengthy  computational  times  required.  In  order  to 
accurately  simulate  a multidither  COAT  servo  and  blooming,  a full  beam 
propagation  calculation  needs  to  be  performed  at  each  computational  incre- 
ment of  the  servo  code.  This  means  roughtly  400  propagation  computations 
during  a single  servo  convergence  cycle.  Since  each  computation  costs 
$2  to  $3,  the  cost  is  prohibitive.  As  a result  of  this  problem,  we  have 
performed  no  simulations  of  multidither  compensation  for  blooming.  One 
way  to  avoid  this  difficulty  is  to  use  a target  irradiance  maximization  scheme 
rather  than  attempting  to  simulate  the  full  multidither  COAT  servo.  We 

plan  on  doing  this  in  the  future,  but  have  not  had  the  opportunity  to  implement 
such  a computation  on  this  contract. 

We  are  beginning  to  accumulate  evidence  that  the  inability  of  a phase 
conjugate  system  to  correct  for  blooming  - and  its  tendency  to  decrease  the 
target  irradiance  with  strong  blooming  - may  be  common  to  all  "return- 
wave"  COAT  systems.  A return-wave  COAT  system  is  one  that  operates 
only  on  the  wavefront  returned  from  the  target  to  derive  the  servo  correc- 
tion signals.  Phase-conjugate  systems,  return-wave  multidither  systems 
and  hybrid  "TRIM-COAT"  systems  are  all  return-wave  systems.  In 
tact,  of  the  presently  known  adaptive  systems,  only  the  multidither, 
outgoing-wave  COAT  system  appears  to  be  a true  maximization  system. 
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We  have  observed,  as  have  others,  that 
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Our  simulation  work  has  also  indicated  that  referencing  off  cf 
extended  target  regions  in  the  presence  of  blooming  will  give  incorrect 
correction  results.  The  problem  is  the  thermal  blooming  equivalent  of 
the  finite  isoplanatic  patch  size  encountered  in  trying  to  image  extended 
targets  through  turbulence.  If  the  target  reference  region  exceeds  an  atmo- 
spheric  coherence  diameter,  a COAT  system  will  not  be  able  to  obtain  the 
correct  phase  errors  for  the  focused  propagation  path.  We  have  observed 
during  our  computer  simulation  work  with  phase  conjugate  correction  for 
blooming  that  this  problem  occurs  even  with  glints  less  than  one -half  the 
diameter  of  the  diffraction -limited  transmitted  beam  on  the  target.  We  plan 
to  study  this  problem  further  on  other  related  programs. 


'The  distance  between  two  points  for  which  two  rays  from  a point  source  see 
the  same  optical  path  length. 
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III. 


PLANS  FOR  THE  NEXT  QUARTER 

The  final  contract  quarter  wil)  be  devoted  p-imarily  to  experimental 
measurements.  We  will  use  the  artificial  turbulence  generators  and  the 
flowing  gas  cell  to  look  at  simultaneous  blooming  and  turbulence  compensation 
for  moving,  multiple -glint  targets.  The  performance  of  the  18  channels  of 
phase  control  will  be  compared  to  that  of  just  the  tracking  and  focus  controls 
and  also  the  performance  of  all  21  channels  operating  simultaneously. 

For  the  computer  simulation  studies,  we  will  begin  looking  for  tech- 
niques that  will  allow  us  to  model  the  operation  of  a multidither  COAT  system 
when  thermal  blooming  distortions  are  present.  The  speckle -noise  simula- 
tion will  be  expanded  to  include  additive  noise  effects  and  we  will  modify  the 
code  so  that  it  can  be  used  on  the  remote  batch-processing  facility  at  HRL 
that  ties  to  the  CDC  7600  machine  at  the  Lawrence  Berkeley  Laboratories. 

We  also  plan  to  run  some  simulation  cases  for  comparison  with  the  acousto- 
optic modulation  experiments  discussed  in  this  report.  The  analysis  technique 
discussed  in  Section  II-B  will  be  used  to  explain  the  computer  simulation  data 
on  speckle  modulation  effects. 
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